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Abstract. Using the IRAM 30 m, SEST 15 m, and Nangay 
radiotelescopes, we have gathered the 1 mm continuum emis- 
sion, the intensities of the J=l-0 line of the CO molecule and 
of the atomic hydrogen line at 21 cm for two samples of IRAS 
galaxies. The southern sample was selected from the IRAS Cat- 
alogue and is complete at the limiting flux of 2 Jy at 60 fim; of 
the 10 northern objects 7 belong to the Smith et al.'s complete 
sample (1987) and 3 are isolated objects. Using these data, 
we have estimated the atomic hydrogen masses from the 21 
cm emission, the molecular gas masses from the CO (1-0) line 
brightness, and the dust and gas masses from the mm con- 
tinuum emission using two "extreme" dust models. The main 
conclusions of this work for far-infrared selected galaxies can 
be summarized in the following points: 

(1) the median value of Mh 2 /Mhi is 0.5, meaning that the 
atomic phase dominates in these galaxies. The fraction of gas 
in molecular form increases with increasing FIR luminosity but 
does not show any obvious trend with other galaxy properties, 
in particular with the FIR surface brightness. 

(2) the H2 surface density derived from CO (1-0) emission is 
better correlated with the cold dust surface density than the 
HI surface density, but the correlation of HI with dust is not 
negligible (we found a correlation coefficient of 0.5, while the 
correlation coefficient with <th 2 is 0.70). Thus, globally in these 
galaxies, the cold dust emission is likely associated with both 
the molecular and atomic phases. Indeed, the dust surface den- 
sity is also correlated with the total gas surface density. 

(3) the FIR surface brightness increases as the third power 
of the S(60 fim) /S(100 fim) ratio. It shows a tight correlation 
with both the H2 and dust surface densities and a weaker one 
with the HI surface density. This suggests that a large part 
of the far-infrared emission of these galaxies originates in the 
molecular medium. 

(4) the gas-to-dust ratio, (Mh 2 + Mm)/Md ranges between 
100 and 1000 and its average value is 230, close to the Galactic 
value. There is indeed a clear trend: this ratio decreases as the 
FIR surface density increases. This result can be explained in 
the framework of an enhancement of metallicity in galaxy discs 
having a higher star formation rate. 
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1. The gas content of galaxies 

The atomic hydrogen phase of the interstellar medium can be 
traced directly from the ground-state emission of atomic hy- 
drogen at 21 cm. But as for its molecular phase, one has to 
rely on more indirect methods. The first one is the observa- 
tion of the CO rotational line (J=l-0), which is supposed to 
trace the column density of molecular hydrogen, and the sec- 
ond is the continuum emission due to dust re-radiation of star 
light in the far-infrared (FIR), submillimeter and millimeter 
wavelength range. The FIR emission is in fact associated with 
hydrogen in atomic, ionized and molecular form. The amount 
of dust is best estimated by combining FIR data with submm 
and mm observations. Indeed, FIR measurements alone, such 
as those of the IRAS satellite, are not sufficient to estimate 
the total amount of dust since they are not sensitive to dust 
colder than 20 K, while an important fraction of the dust emis- 
sion is expected to lie at submillimetre /millimetre wavelengths. 
Moreover, the evaluation of the physical properties of the dust 
grains are less affected by observational uncertainties if obser- 
vations are carried out in the submm/mm range, where the 
relationship between the flux and the dust temperature and 
mass is linear . 

All these procedures provide an estimation of the gas mass 
present in the associated medium, but it must be noted that, 
while the column density of the emission line at 21 cm gives a 
direct evaluation of the atomic hydrogen mass, the conversion 
factors between the dust and CO emissions on the one hand, 
and the gas mass on the other hand, depend on the physi- 
cal conditions in the galaxy discs, conditions which are not 
very well known. The comparison between the estimates of gas 
masses from the CO emission and the mm continuum emis- 
sion should reflect the large uncertainties involved in both ap- 
proaches, but, since they are based on completely independent 
methods, this comparison would be of great help in assessing 
their reliability, calibration and, therefore, in discriminating 
among models of the interstellar medium. 

Many papers (e.g., Young et al. 1986; Mirabel fc Sanders 
1989; Young et al. 1989; Devereux fc Young, 1990; Tinney et 



Send offprint requests to: P. Andreani 



(*) based on observations collected at the European South- 
ern Observatory, La Silla, Chile 



ai. i»»u; oanaers ei ai. n.aaiora ei ai. uownes, 

Solomon fc Radford, 1993) have already addressed the question 
of the relationships between the CO, HI and FIR emissions 
in IRAS galaxies. Most of them, however, deal with extreme 
sources, such as luminous IRAS galaxies and/or active galaxies. 

The basic results of previous works can be summarized in 
the following points. 

(1) FIR and CO luminosities are strongly correlated. Their ra- 
tio increases with the FIR luminosity and the warm dust tem- 
perature (given by the 60 /zm-100 fim colour 560/5100 )• It has 
been suggested (Tinney et al., 1990) that this ratio is a func- 
tion of the interaction stage: it ranges between 20 for isolated 
objects and ~ 300 for mergers. 

(2) the FIR luminosity correlates better with the molecular 
mass, as derived from the CO emission, than with the HI mass 
or the total HI + H2 mass, indicating that warm dust is pri- 
marily in molecular clouds. 

(3) the mean value between the H2 mass and the warm dust 
mass derived from the far infrared data, (MH 2 /Afwd) is of the 
order of 500, that between the total gas mass (H2 + HI) and 
the warm dust mass, (Af gas /Af w d} is roughly 1000, much larger 
than the galactic values of roughly 160 (see e.g. Sodroski et 
al., 1994), suggesting that the bulk of the dust could lie at 
temperatures of the order of 20 K. 

(4) for Mb_ 2 ranging between 10 8 and 5 10 10 Mq, the ratio 
Mh 2 /Mhi has been claimed to be close to 1 (Young fc Knezek 
1989), but other samples yield much lower values, about 0.1 
(e.g. Sage 1993). 

(5) The higher the FIR luminosity, the larger is the ratio 

Mh 2 /M H i. 

Note that most of the above analysis have been made by 
comparing masses and luminosities of the galaxies. However, it 
is well known that large galaxies are also more luminous at all 
wavelengths, so that in order to derive the true relationships 
between the observed quantities, one has to get rid of the size 
effect by scaling the data to the total mass or the surface of 
the galaxy. 

At long wavelengths, in the submm/mm range, very few ob- 
servations are available in the literature. A comparison between 
the continuum emission at 1.25 mm and that of the CO(l- 
0) line for a sample of IRAS Markarian galaxies has been per- 
formed by Kriigel et al. (1990), Chini et al. (1992), and Chini 
fc Kriigel (1993). These authors discuss their data in terms of 
star formation efficiency given by the ratio between the FIR 
luminosity and the gas mass derived from the cold dust column 
density, Xfir / -Mgas . They find that the mean value of this ratio , 
(ipiR/Afg as ), is 123 ± 56. They also find a strong correlation 
between ipiR/Af gas and the cold dust temperature. The edge- 
on galaxy NGC 891 has been mapped at 1.3 mm by Guelin et 
al. (1993) : the spatial distribution of the continuum emission 
matches very well the distribution of CO (1-0) emission and 
the gas masses derived from the dust and CO emission agree 
quite well. 

These results suggest that by combining FIR and mm mea- 
surements of the dust emission, one should be able to measure 
the neutral gas content of galaxies. The aim of the present work 
is to compare gas masses computed by this method with those 
using HI and CO data, and to use this comparison to derive 
some of the physical properties of the interstellar medium in 
normal galaxies. 

Section 2 deals with the three sets of observations, Sect. 
3 reports the procedures to compute the gas masses from the 



ODservea nuxes. ine statistical analysis is aescriDea in oec- 
t. 4, while the conclusions drawn from the present work are 
summarized in Sect. 5. 



2. The observations 

We have observed two samples of galaxies : a northern one, 
for which the millimeter observations have been done with the 
IRAM 30 m telescope in Spain, and a southern one, for which 
the observations have been done with the SEST 15 m telescope 
in Chile. Galaxies observed with SEST belong to a complete 
sample of objects selected from the IRAS Point Source Cata- 
logue, for which optical photometry and spectroscopy are avail- 
able. The sample includes 29 galaxies with 60 fim fluxes above 
2 Jy in the sky region delimited by the equatorial coordinates 
21 h < a < 5 h and -22.5° < 6 < -26.5°. 

Due to the small beam size and in order to cover as far as 
possible the entire galaxy and to reduce the effects of aperture 
corrections, observations at IRAM were limited to 7 northern 
galaxies (8 h < a < \1 h and 23.5° < 8 < 32.5°) from the 
Smith et al. sample (1987) plus 3 isolated objects, with optical 
angular dimensions smaller than or comparable to the beam 
width. This may have introduced a bias relating to the source 
compactness in the results obtained for this sample. 

Table 1 reports name, coordinates, distances, 25th magni- 
tude diameter, B magnitudes, velocities, morphological type, I- 
RAS fluxes and FIR luminosities for the observed galaxies. The 
data have been taken from the Lyon-Meudon Extragalactic 
Database (LEDA), the NASA Extragalactic Database (NED), 
and IRAS PSC except for some objects of the northern sample, 
which have been taken from Smith et al.'s paper (1987). The 
FIR luminosities have been computed by adding the 1.3mm 
contribution to the luminosity in the wavelength range 40 - 
120 fim, calculated as in Lonsdale et al. (1985). 

Throughout this paper, we take H = 75 kms" 1 Mpc" 1 . 
Note that all velocities quoted here are heliocentric, using the 
optical definition v op t = cz = c AA/A. 



Table 1 : the sample 
2.1. The bolometer observations 

The continuum observations at 230 GHz have been taken dur- 
ing March and May 1990 with the IRAM 30 m telescope, and 
during September 1990 and 1991 with the SEST 15 m tele- 
scope. Observations have been already described elsewhere 
(Andreani et al. 1990, Andreani fc Franceschini 1992; Frances- 
chini fc Andreani, 1994), here we briefly outline the observa- 
tional set-up. Both antennas fed a 3 He-cooled bolometer of the 
MPIfR (Kreysa 1990), whose filter set, coupled to the atmo- 
spheric transmission window, provides an effective wavelength 
around 240 GHz (1.25 mm) with a bandwidth of 50 GHz. The 
beam size was 11" and 24" (HPBW) for IRAM and SEST re- 
spectively, while 

the wobbling amplitude and the chopper throw were set 
commensurate to the average optical dimensions of the galaxies 
(i.e. 50" and 70" for the IRAM and SEST, respectively). 

Since most of the objects observed with SEST and IRAM 
have optical sizes comparable to the telescope beam aperture, 



aperture corrections may men turn out to ae signmcant ana de- 
serve particular attention. Moreover, the beam sizes at 1.25 mm 
are twice smaller than that at 2.6 mm and much smaller than 
that of the Nangay antenna (4' x 21'). This means that if the 
source extension at long wavelengths is comparable to the op- 
tical dimensions (see Table 1) an important fraction of the 
1.25 mm flux could have been missed. To bypass this problem, 
scans along the major axis have been performed for the most 
extended sources (actually NGC1187, NGC1385, NGC7252, 
NGC7314). The fluxes have been obtained by adding the con- 
tributions of each observed position. For the remaining sources, 
for which scans are not available, a small correction for beam- 
aperture has been applied. Following Andreani fc Franceschini 
(1992) and Franceschini fc Andreani (1994), we assume that 
the radial distribution of the millimetric light is exponential 
with a scale-length a mm equal to one third of the optical one: 
c*mm = c*o/3. This means that the dust emission is more con- 
centrated than the star-light. Aperture-correction turns out to 
be of the order of 10-30% for a few sources, namely NGC7225, 
E534-G9, NGC578, E478-G6, NGC922 and negligible for the 
remaining ones. 

Nineteen objects have been detected with SEST and for 
another nine an upper limit at 3 a has been given. Of the 
nine IRAM sources, five have been detected and for the re- 
maining five an upper limit at 3<r has been put. The data on 
the 1.25 mm continuum emission have been already published 
elsewhere (Andreani et al. 1990, Andreani fc Franceschini 1992; 
Franceschini fc Andreani, 1994) and the complete observation- 
s will appear together with a wide discussion on them in a 
forthcoming paper (Franceschini fc Andreani, in preparation). 



2.2. CO Observations 

We have observed the CO (1-0) emission of the southern ob- 
jects in our sample in November 1992, using the Swedish-ESO 
SubmiUimeter Telescope (SEST) in La Silla (Chile). The beam- 
size is 43" and the main-beam efficiency is 0.74. We used a 
Schottky receiver in single sideband mode with T rec PS 300 K ; 
the system temperature in the TjJ scale varied between 400 and 
500 K. The backend was an 1440 channel acousto-optical spec- 
trometer with a channel width of 691 kHz, but to improve the 
signal-to-noise ratio the data were smoothed to a resolution of 
22 kms - . Throughout the article, the brightness temperatures 
will be expressed in the main-beam scale, T m b = T^/i/mb. The 
conversion factor from antenna temperatures to flux densities 
is S;,/T m b = 27 JyK - at 115 GHz for this telescope. 

The observing procedure was dual beam-switching at a 6 Hz 
rate with two symmetric reference positions offset by 12' in az- 
imuth. The observational strategy was the following : for each 
galaxy, we made an estimate of the expected CO emission from 
its far-infrared emission, using the relationship Xpir/Mhj = 12 
and a conversion factor from CO emissivities to H2 column 
densities of 2.310 cm - (Kkms - ) - (note that this value 
will be used throughout the paper to derive H2 masses). We 
have then estimated the integration time necessary to reach a 
signal-to-noise ratio of 3, for a linewidth of 300 kms -1 . The 
central position was observed for the time determined in this 
way. For eight galaxies with an optical diameter D25 (as giv- 
en in the RC3) larger than than 110", we have made maps at 
45" sampling along the major and minor axes. The total in- 
tegration time on the off-center positions was the same as for 



me central point, vvim mis procedure, me typical integration 
times ranged between one and two hours. 

Out of the 25 galaxies observed, 19 were detected (the 
detection rate is then 80%). One galaxy (NGC232) was not 
detected by us, but it has been observed with the SEST by 
Mirabel et al. (1990) and detected. The data for NGC 7252 
come from our previous observations with SEST (Dupraz, Ca- 
soli fc Combes, 1990), while those for IC860 come from our 
CO (1-0) survey of the Coma cluster (Casoli et al, in prepa- 
ration). Data are listed in table 2: names are in column 1, 
CO (1-0) line intensities (7co given in Kkms) and velocities 
are given in column 5 and 6, velocity widths in column 7. 

2.3. HI Observations 

Since we wanted to compare dust masses with total gas masses, 
i.e. molecular plus atomic hydrogen, we have obtained HI da- 
ta for those galaxies which have never been observed at 21 cm. 
They were observed with the Nangay radiotelescope (France) in 
1993. The half-power beamwidth of the instrument is 4' (E-W) 
by 21' (N-S). The dual channel receiver has a system temper- 
ature around 50 K. The backend was a 1024 channel autocor- 
relator. Total power mode was used with a reference position 
30' away from the source, and the total integration time var- 
ied from one to four hours. For each polarization, the antenna 
temperatures have been converted to flux densities using a con- 
version ratio which depends on the polarization of the receiver 
and on the declination of the source ; it is typically 1.1 JyK -1 
for the H polarization and 0.8 for the V polarization. The fi- 
nal spectra were obtained by addition of the two polarization 
modes, smoothing to 10.55 kms -1 and subtraction of a linear 
baseline. In total, 28 galaxies were observed and 23 detected 
(out of the detected ones, there is E478-G6 for which there 
are clearly two galaxies in the Nangay beam; we have chosen 
the line whose central velocity is closest to the CO and optical 
one). The HI data for other galaxies has been extracted from 
various data bases (the NASA Extragalactic Database, NED 
and the Lyon-Meudon Extragalactic Database LED A). All the 
HI data are listed in table 2:21 cm line fluxes (in Jy) and ve- 
locities (in kms) are given in column 2 and 3, velocity widths 
(in kms) in column 4. 



Table 2 : CO and HI data 



3. Determination of dust and gas masses 

3.1. Gas masses from dust models 

As the thermal dust emission at mm and submm wavelengths 
is optically thin, at these wavelengths the entire line of sight is 
sampled and the total dust column density can be determined. 
The relationship between the dust mass Md and temperature 
Td (it is possible to approximate the temperature distribution 
of the grains with a single temperature), is then linear: 

1.631 10 8 A 2 (cm)D 2 (Mpc)f A (Jy) 

AZ d (M ) - 2 _ n7 , /^x 

XA(cm 2 g 1 )T d (K) 

where %a is the grain opacity and D is the source distance. 



ine main uncertainly in mis relation lies in me meorencai 
assumptions, since different dust models predict significantly 
different values (up to a factor of 10) for the longwavelength 
grain opacity. We have estimated the amount of dust and gas 
by using two separate (and "extreme") models for the inter- 
stellar dust (Draine fc Lee 1984, hereafter DL84; Mathis fc 
Whiffen 1989, hereafter MW89), which assume different com- 
position and grain dimensions. Dust and Gas masses have been 
estimated from both dust models using for MW89 a grain opac- 
ity at 1.2mm of 1.74 cm per g of dust and a gas-to-dust ratio of 
366, while for DL84 the grain opacity is 0.225 cm 2 per g of dust 
and the gas-to-dust ratio 159. The dust (colour) temperature, 
given by the S(100/zm)/S(1.25mm) ratio, ranges between 15- 
25 K for DL84's model and it is 30% lower than that resulting 
from MW89's model. 

In the following, gas/dust masses and gas/dust surface den- 
sities calculated with MW89's model are denoted as M* 1 /M^ 1 
and (jg 1 ' /cjj"' while those calculated with DL84's model are de- 
noted as Mf ] /Mf ] and a [ g ] /a^ ] . 

These two dust models have been chosen because their predict- 
ed grain opacities (per g of dust) at long wavelengths differ by a 
factor of 7. This difference is taken as a conservative confidence 
range for the actual dust mass values. 

Table 3 lists the dust and gas masses from MW89 model 
(column 2 and 4) and DL84 model (column 3 and 5) for the 
observed objects. For those not having detections a 3<r upper 
limit has been given. 

3.2. H2 masses 

The H2 column densities were obtained from the CO (1-0) line 
areas using a conversion factor of 2.310 cm - (Kkms - ) - 
(Strong et al. 1988). For the SEST 43" beam, this converts 
into : 

Mh 2 (M ) = 1.25 10 5 7 co (Kkms" 1 )(D/Mpc) 2 

in the beam, where D is the distance to the galaxy in Mpc 
and 7co is the integrated CO line intensity. For those readers 
who prefer to use CO luminosities, the relationship between 
the H2 masses given in Table 3 and the CO luminosities is 
Mh 2 (M ) = 3.67ico(Kkms _1 pc 2 ) . Note that this value a- 
grees with that used by Radford et al. (1991) but it is smaller 
than that adopted by other authors: Sanders et al. (1991) uses 
a value higher by a factor of 1.3, while Young et al. (1986) a 
value higher by a factor of 1.6. The value we adopted for the 
conversion factor has been derived by Strong et al. (1988) from 
the analysis of the galactic 7 rays emission measured by the 
COSB mission. 

3.3. HI masses 

We have computed the HI column density in the beam using 
the relationship : 

Mhi(M ) = 2.36 10 5 S V AV(D/1 Mpc) 2 

where S„AV is the line flux in Jykms -1 . Table 3 gathers the 
derived HI and H2 masses (columns 6 and 7). 

Table 3 : dust and gas masses 
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4-1. The Survival Analysis for univariate censored data 

In order to make a full use of the information contained in the 
upper limits to the 1.25 mm, CO and 21cm fluxes, we have 
used the technique of survival analysis. The Kaplan-Meier al- 
gorithm (see e.g., Feigelson fc Nelson 1985; Schmitt 1985) has 
been used to estimate the parent distributions and the average 
values of the gas masses for the three different sets of obser- 
vations and of the ratios between the gas mass found from 
the dust column density and the molecular gas mass given by 
the CO luminosity, /Mu 2 and M^ 2I /Mh 2 , for the two dust 
models, as well as the ratio between the molecular and atom- 
ic gas masses, Mh 2 /Mhi- We have also computed the average 
values of the logarithm of the gas masses, which are more repre- 
sentative of the most probable value because of the large range 
and asymmetry of the observed quantities. We find as mean 
values (we also give the estimated error of the mean): 

(Mf 1 ) = (1.9±O.5)1O 9 M , (log(M[ 1] )) = log(0.9 1O 9 M )±O.1, 

(M| 2] ) = (7.7±2.O)1O 9 M , (log(Mf ] )} = log(3.1 1O 9 M )±O.1, 

(Mh 2 ) = (3.1±1.O)1O 9 M , (log(M H2 )> =log(1.3 1O 9 M )±O.1, 

(Mm) = (3.7±O.5)1O 9 M , (log(Afm)) = log(2.30 1O 9 M )±O.O7. 
The median values are : 

med(Af| 1] ) = 1.O1O 9 M , 

med(Af| 2] ) = 4.1 10 9 M , 
med(MH 2 ) = 1.2 10 9 M , 
med(Mm) = 2.5 10 9 M , 

The meaningful variables are actually the gas surface densi- 
ties, i.e. the ratios of the gas mass to the optical surface of 
the galaxy given by D| 5 /4 where D25 is the 25-th magnitude 
diameter given in the databases. Note that these are hybrid 
surface densities since they mix optical and gas parameter- 
s. Ideally, to compute the HI surface density for example, one 
would like to use the HI characteristic diameter which is known 
to be larger than the optical one in most galaxies, but this 
quantity is generally unknown. As for the CO and dust e- 
missions their scale lengths are generally smaller than that 
of the visible light, but also unknown for the sample object- 
s. However, we assume that the ratio between the dust/CO 
scale-length and the optical one remains constant for this kind 
of objects. Indeed, Haynes and Giovanelli (1984) have shown 
that the HI surface density is roughly constant at the value 
(log(<JHi)} = log(6.5 lO 6 )M /kpc 2 ) in spiral galaxies and does 
not present any strong trend with the galaxy luminosity. Such 
surface densities are thus useful parameters since they are in- 
dependent of the assumed distance of the galaxy, on the one 
hand, and show weak dependence on the galaxy size and mor- 
phological type, on the other hand. 
The values we find are: 
(log(4 1] )) = log(2.1 10 6 M /kpc 2 ) ± 0.1, 

(log^ 21 )} = log(8.4 10 6 M /kpc 2 ) ± 0.1, 

(log(<7H 2 )> = log(2.5 10 6 M /kpc 2 ) ± 0.1, 

(log((7Hi)) = log(5.3 10 6 M /kpc 2 ) ± 0.1, 
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For these five variables, the median values are similar to 
the averages; we found respectively : log(1.710 6 ), log(7.810 6 ), 
log(2.5 10 6 ), log(5.2 10 6 ) and log(7.4 10 6 ). 

The average HI surface density for these IRAS selected 
galaxies is thus close to the canonical value for isolated galax- 
ies. As for log(<JH 2 ), the analysis of the FCRAO extragalactic 
survey performed by Gerin fc Casoli (1994) yields a value of 
(log(<JH 2 )} = log(2.2 10 MQ/kpc ). This value is very close to 
the one derived above, which probably reflects the fact that 
the FCRAO sample contains mainly galaxies that are bright 
either in the infrared or in the optical. On the other hand, the 
distance-limited sample studied by Sage (1993) shows a much 
lower average H 2 surface density of log(0.710 6 M /kpc 2 ). The 
total gas (H2 + HI) surface density computed from CO and HI 
shows a remarkable constancy from one sample to the other 
being log(8.3 10 6 M /kpc 2 ) for both the FCRAO and Sage's 
samples (Gerin fc Casoli 1994). The IRAS selected galaxies 
show the same average gas surface density. 

Fig.la ; 6 shows the cumulative distributions for the ratios 
between gas masses determined from dust models and from CO 
and HI measurements: log(Afg 1 '/Af gas ), log(Afg 2 '/Af gas ), where 
Mgas is the total gas mass derived from the CO and HI data : 
Mgas = Afai + Afa 2 • Note that each value of the distribution 
represents the total probability for the population to have a ra- 
tio lower than the corresponding one on the x-axis. The arrows 
indicate the average values: 

(log(M g 1] /M gas )> = log(0.28), 

(log(M g 2] /M gas )> = log(1.33). 
The total gas masses and surface densities are slightly lower 
than that inferred from DL84 model and a factor of 3 higher 
than that found from MW89 model. This result is not com- 
pletely unsatisfactory if one considers that the procedures of 
estimating the gas masses are completely independent and for 
both the associated uncertainties are quite large. 
The corresponding gas-to-dust ratios are : ((M gaiS /M^)) ~ 
1300, and ((M gas /M| 2] )> ~ 120. 

Panels c,d,e,f of the same figure report the cumulative 

distribution of the ratio between gas masses: log(Af g 1 '/MH 2 ), 

log(M g 2] /M H2 ), log(Af Ha /Afm) and log(M H2 /(M H2 + M m )), 
with fairly well-defined average values: 

(log(M g 1] /MH 2 )> ~log(0.7)±0.1, 

(log(M g 2l /M H2 )> ~log(3.2)±0.1, 
(log(M H2 /M H i)) ~ log(0.53) ± 0.10, 
(log(M H2 /M gas )> ~ log(0.30) ± 0.06. 

4-2. Regressions with censored data 

Following Schmitt (1985), we used his procedure to perform 
linear regression when both the independent and dependent 
variables have upper limits (see also Isobe, Feigelson and Nel- 
son, 1986). 

Let us call the variables (X,Y) and the total number of 
points iVtot. The (X,Y) plane has been divided into (N,M) 
bins (JV the number of bins for the dependent variable and 
M that for the independent variable) and to each bin (xi,yj) 



Fig. 1. Cumulative distributions of the following ratios between gas 
masses: (a) log(M g 1] /M gas ), (b) log(M g 2] /M gas ), (c) log(M g 1] /M H2 ), 

(d) log(M g 2l /M H2 ), (e) log(M H2 /M m ), (f) log(M H2 /M gas ). Arrows 
indicate the average values reported in the text (§4.1). [1] corre- 
sponds to predictions of the Mathis and Whiffen's dust model; [2] 
to those of the Draine and Lee's model. 

a two-dimensional probability density function similar to the 
univariate Kaplan-Meier function is assigned, f(xi,yj). Once 
the bins have been ordered from the largest to the smallest one 
and in each bin the four different types of data (detected points 
D(x,y), points with a censored independent variable U x (x,y), 
points with a censored dependent variable U y (x,y) and points 
with both censored variables U xy (x, y)) have been counted, the 
two-dimensional probability density function is estimated by 
iteratively solving the maximum likelihood equation given by 
Schmitt: 

El-i t/ x(a:i,2/)/(a;,2/) 
N tot f(x,y) = D(x,y) + — + 

Once the probability density has been estimated, correlation 
coefficients and linear fits are obtained by taking the various 
moments of the probability density. It is not straightforward 
to compute the confidence range of the regression values. Er- 
rors are computed by bootstrapping the data: many different 
samples, each containing Ntot data as the real one, are created 
by randomly redistributing the data from the real catalogue, 
i.e. the sequence number of the data is randomized, this means 
that in some cases one datum can be considered more than 



once, wmie in some omer cases 11 can ae aiscaraea. ine proD- 
ability density function is computed for each new sample and 
the estimated uncertainties is taken as the standard deviation 
from the mean value of the function, value obtained by aver- 
aging the functions over all the created samples. 

Results of the regressions are shown in fig. 2^6. The solid 
lines are the computed regression lines. In each panel are are 
given the values of the correlation coefficient, p, and of the 
slope with its uncertainty. 



5. Discussion : comparison between the various 
tracers of the gas content 

5.1. Correlation dust - gas 

The two adopted dust models assume quite 'extreme' dust 
properties at long wavelengths and therefore their predictions 
differ by a large factor. We think that this difference can be 
taken as a very conservative confidence range for the predict- 
ed dust properties. We use, therefore, in the following as dust 
masses the (linear) average values of the ones given by the two 
models. Fig. 2 shows the behaviour of the H2 and HI surface 
densities and of the total (H2+HI) one as a function of the dust 
surface density. From an inspection of the figure the following 
conclusions can be drawn: 



Fig. 2. Regression between the logarithms of the following couples 
of variables: (a) [cth 2 , <r d ], (b) [a m , <7 d ], (c) [<t H i+h 2 , <r&], ( d ) [<m 2 , 
<thi]- The solid-line corresponds to the computed regression line. The 
related values of the correlation coefficient and the regression slope 
are reported in the top left corner of each panel. 



(1) the HI surface density, the H2 surface density derived from 
the CO emission, as well as with the total gas surface density 
CTHI+H2 are correlated with the dust surface density computed 
from the mm-continuum emission. The best correlation is with 
<th 2 , the weakest one is with <thi, but none of these correlations 
is excellent. Because the slopes are lower than 1, (they range 
between 0.6 and 0.8), the gas surface densities are not pro- 
portional to the dust surface densities. This suggests that the 
gas-to-dust ratio varies from one object to another. From these 
data we can conclude that the cold dust in these IRAS-selected 



is contained Dom in moiecuiar ciouas ana in me atomic medi- 
um, but more closely associated with the molecular gas. 

The continuum emission can then be used as a tracer of the 
total hydrogen column density (atomic and molecular) insid- 
e the optical radius, independently of the CO (1-0) emission, 
but the uncertainties are still large. Note that the observation 
time required to detect the galaxies at 1 mm is similar to that 
necessary to reach the detection of the CO emission. However, 
the detection rate is slightly better in CO than in the mm- 
continuum. 

(2) Panel (d) of figure 2 shows the behaviour of <jh 2 versus <thi. 
The correlation is weak (p = 0.5) and shows a large scatter. 
On average the HI content in these galaxies is twice that of H2 
(see §4.1). 



5.2. H2 gas fraction 

Fig. 3 shows the correlation between the H2 gas fraction and 
the FIR luminosity, (the correlation coefficient is 0.75). This 
is clearly linked to the fact that Mh 2 and £fir are correlated, 
while neither Mh 2 and Mm, nor Mm and Xfir, are strongly 
correlated (see Fig. 2d and Fig. 4b). The behaviour of the log- 
arithm of Mh 2 /(Mh 2 + Mm) against that of the FIR colour, 
5(60)/5(100) and of the FIR surface brightness, ctfir, indi- 
cates a trend of increasing H2 gas fraction with increasing ctfir 
and 5(60)/5(100) (correlation coefficients 0.51 and 0.53). We 
do not see any strong trend of the molecular gas fraction with 
the S(100)/S(1250) colour, thus with the cold dust tempera- 
ture. This could be interpreted in view of the results of figure 2: 
cold dust is associated with both molecular and diffuse media. 



Fig. 3. Relationship between the logarithmic values of the quanti- 
ties: [Mn 2 /(Mh 2 + .Mm) and Lfir]- The solid-line corresponds to 
the locus where objects having Mn 2 = Mm lie. 



It is interesting to note that the solid line in Fig. 3 cor- 
responds to the case Mm = Mh 2 . Most objects lie below this 
line, having Mh 2 < Mm. The figure clearly points out that the 
H2-rich objects are all bright in the far-infrared, having FIR 
luminosities larger than 310 10 Lq. Several mechanisms have 
been advocated to explain a possible depletion of HI gas in IR- 
bright galaxies (ionization of atomic hydrogen because of the 
intense UV field and/or induced by supernova winds) but the 
more likely process is the conversion of atomic to molecular gas 
: since the largest part of the gas mass is in the atomic phase, 
the conversion of a small fraction of this atomic gas into molec- 
ular gas leads to a large variation of the fraction of molecular 



gas dui leaves me mass 01 atomic gas nearly uncnangea. ±ms 
mechanism could be enhanced in merger systems (Mirabel fc 
Sanders, 1989). 

However, it is possible that in some objects a fraction of 
the H2 gas is not traced by the CO emission. This is true for 
example in low-metallicity environments (for instance in the 
outer galaxy or in the Magellanic Clouds), where CO forms in 
the deepest regions of molecular clouds with respect to those 
where H2 does, or in very cold gas where the rotational exci- 
tation of CO is weak (Lequeux, Allen and Guilloteau, 1993; 
Pfenninger, Combes and Martinet, 1994). We do not expect 
however our galaxies to be very metal-poor since they must 
contain a reasonable amount of dust in order to be detected 
in the FIR and submillimeter range. Young fc Knezek (1989) 
claim that the Mh 2 /Mhi ratio is higher in early-type spirals. 
At first glance this seems not to be the case for our sample, 
however we lack any statistical significance for this test. 



5.3. Comparison between arm and ad, cth 2) cthi, 60/100 

Further investigation of the relationship between FIR luminos- 
ity and dust and gas masses can be pursued by comparing the 
FIR surface density with the H2 , HI and dust surface densities 
and also with the warm dust colour. 



Fig. 4. Relationship between the logarithmic values of: (a) [<xfir 1 
<m 2 ], (b) [-TFIR , <thi], (c) [-TFIR , <T d ], (d) [<X F IR , 5(60)/5(100)] 
The solid-line corresponds to the computed regression line. The re- 
lated values of the correlation coefficient and the regression slope are 
reported in the top left corner of each panel. 



Figure 4a shows a plot of log(crpiR) against log((JH 2 ), which 
confirms that the well-known correlation between FIR and CO 
luminosities also holds for the surface densities. This tight rela- 
tionship (p = 0.87) is more remarkable since both variables are 
(to the first order) independent of the distance and of the sizes 
of the objects and are thus really intrinsic properties. Note that 
contrary to most of the correlations studied in this paper, this 
relationship is linear within the errorbars: 

log((7FIR) = (1-14 ± 0.15) log((7H 2 ) + (0.4 ± 0.9) 



inere is a weaxer correlation 01 ctfir wnn cthi mg.^Dj. ims 
means that the FIR emission associated to atomic clouds is 
weak compared to other sources in IRAS galaxies. There is a 
good linear relationship between ctfir and <Jd (fig. 4c), thus 
confirming that most of the millimetric emission is thermal. 

It is interesting to investigate the relationship between 
log(crpiR) and the logarithm of the ratio between the 60 and 
100 (im flux densities, log(560/5100), plotted in fig.4d: most 
of the objects follow the relation: ctfir cc (560/5100) (dashed 
line in fig.4d), with correlation coefficient of 0.7. The increas- 
ing of both dust temperature and FIR surface brightness as 
the UV radiation field increases could give rise to the strong 
variation of ctfir with the FIR colour seen in figure 4d (see 
also Sauvage and Thuan, 1994). 



5.4- The Lyir/Mb_ 2 ratio in our sample 

We have compared the Xfir /Mh 2 ratio with Xfir, Mh 2 
and the ratio between the 60 and 100 am flux densities, 
S(60)/S(100). Fig. 5 shows that our sample does show nei- 
ther a dependence of the Xfir /Mh 2 ratio on Xfir nor on 
Mh 2 . The median value of log(i F iR/M H2 ) is log(9.44) ± 0.14 
, typical of IRAS galaxies (e.g. Sanders et al., 1986; Sanders 
et al., 1991). Only two nearby objects (NGC1187,NGC1255) 
have Lyir/Mb_ 2 ~ 3Lq/Mq, while most lie in the range 8-18 
Lq/Mq. The galaxy IC860 is another exception, which has 
Lyir/Mb_ 2 ~ IOOXq/Mq, but it is a compact object and be- 
longs more likely to the class of peculiar galaxies. Indeed it ex- 
hibits both OH megamaser and H2CO maser emission (Baan 
et al 1993). 

Young et al. (1986) have found a good correlation be- 
tween the ratio of FIR and CO luminosities, and the warm 
dust temperature. Figure 5c shows the comparison between 
log(iFiR/MH 2 ) against the logarithm of the ratio of the 60 
and 100 am flux densities. There is indeed a trend of increas- 
ing Lyir/Mb_ 2 as S(60)/S(100) increases (with a correlation 
coefficient of 0.58) which could be interpreted as an evidence 
for FIR emission being mainly thermal. 

Panel (d) of figure 5 shows a plot of the ratio between FIR 
luminosities and total gas masses Mh 2 + A/hi, £FHi,/M 3as , as 
a function of the FIR luminosity. A correlation between these 
two quantities is seen (p = 0.68). The reason is that, while 
Mh 2 and Xfir are roughly proportional, Mm does not depend 
strongly on the FIR luminosity and thus the total mass of 
neutral gas M gas does not increase linearly with Xfir but with 
an exponent smaller than unity. 

Chini et al. (1992) found for a sample of active galaxies 
an average value of the (iFiR/M gas ) ratio of 123 ± 56, where 
M gas is the gas mass found from dust emission and Xfir is the 
total FIR luminosity from 12 to 1300 am. Our sample has an 
average value of log(iFHi,/M gas ) of log(5.4)±0.46 (the median 
value is 4), where M gas comes from Mh 2 + Mm. If we use 
the gas mass inferred from dust emission, we still find a value 
between log(3) and log(ll), one order of magnitude lower than 
that of Chini et al.'s. This difference cannot be accounted for 
by the different ways of computing the far-infrared emissivities 
, which could amount to a factor of 2, and it could be ascribed 
to properties related to galaxy activity: our sample consists 
mainly of normal galaxies, while that of Chini et al. of active 
galaxies. 



Fig. 5. Comparison of the logarithm of the quantities: (a) £fir/^H 2 
andL FIR , (b) L FIR /M H2 andM H2 , (c) L FIR /M H2 and 5(60)/S'(100)| 
(d) L FIR /(M H2 + M m ) and L FIR . 

5.5. Gas-to-dust ratio 

The gas-to-dust mass ratio has been evaluated for our sample 
using the total gas mass Mh 2 + -Mm and the dust masses es- 
timated from the average value given by the two dust models. 
Figure 6 shows the relationship between the gas-to-dust ratio 
and the FIR surface density. There is a tendency for this ratio 
to decrease as <t F i R increases. The interesting point to note is 
the wide range of the observed values for the gas-to-dust mass 
ratio. For most objects this value is around 300 for MW89 
and 150 for DL84, which are indeed close to the solar neigh- 
borhood values appropriate to these models (see Sect. 3.1). 
The gas-to-dust ratio is very large (>1000) for NGC1187 and 
NGC578, which are large objects (about 5 arcmin in diameter) 
with low Lyir /.Mh 2 ratios and cold FIR colors. The most like- 
ly explanation is that we could have missed some of the dust 
emission or underestimated the correction to be made (sect. 
2.1). On the other hand, it could be argued that these two 
objects are truly metal poor but this is unlikely since we have 
detected CO emission in both of them. Some galaxies appear 
to be very dusty : the gas-to-dust ratio is below or around 100 
for IC 1553, 2101-23, E484-G36 and IC860. Of these galaxies 
only IC860 shows clearly starburst signspots: warm infrared 
colors and high far-infrared luminosity compared to their opti- 
cal ones. The other three objects have indeed FIR luminosity 
larger than the optical one. If this property alone is related to 
star-formation activity, it could be argued that, as other obser- 
vations indicate, this latter increases the metallicity and there 
is indeed a proportionality between the dust-to-gas ratio and 
heavy element abundance (e.g. Issa, McLaren fc Wolfendale, 
1990; Sodroski et al., 1994) This ratio could be used therefore 
as a metallicity indicator. 

It must be noted that on the basis of the 60 and 100am 
fluxes only, it is not possible to infer the dust temperature 
and mass (see Draine, 1990 for a discussion on this point) and 
nothing can be said about the overall gas/dust ratio. The mil- 
limetric points undoubtly show that the 100 am fluxes trace 
better the cold component to which the mm data belong, while 



Fig. 6. Comparison between gas/dust ratio andlog(<r F j R ). The solid 
line corresponds to the computed regression between the two quan- 
tities. In the upper right corner the correlation coefficient and the 
slope with its uncertainty are reported. 

the 60 am fluxes, mainly in starburst galaxies, are a signature 
of the warm component (Chini fc Kriigel 1993, Franceschini 
fc Andreani, 1994). Moreover, the dust masses inferred on the 
basis of the 60/100 colour only (warm dust mass) are on the 
average one order of magnitude lower than the dust masses 
(warm and cold dust) found with the entire FIR/ mm data. In 
addition, temperature fluctuations of the smallest dust parti- 
cles may still be important at 60 am (Desert, Boulanger and 
Puget 1990 ; Sodroski et al 1994). 

5.6. Possible biases 

Our sample of galaxies has been selected on a FIR selection 
criterion, as was the Chini et al.s (1992) one. It is clear that 
galaxies which have a strong far-infrared emission must have 
both a large dust abundance and active star formation. It is 
known that this biases the sample towards CO-rich galaxies. 
The same kind of investigations should be pursued also for an 
optically complete sample of galaxies. Submillimeter observa- 
tions at shorter wavelengths would be of great help. Finally, 
the use of bolometer arrays should help to determine the extent 
of the dust emission. 



6. Conclusions 

The main conclusions of this work are summarized in the fol- 
lowing: 

(1) the median value of (Mh 2 /Mhi) is 0.5. Contrary to the 
what has been observed in the FCRAO extragalactic sample 
(Young fc Knezek 1989), the atomic phase dominates in these 
galaxies. Part of the difference is due to a different conversion 
factor of the CO (1-0) emissivities to the H2 column densities. 
The fraction of gas in molecular form increases with increasing 
FIR luminosity but does not show any strong trend with other 
galaxy properties, in particular with the FIR surface bright- 
ness. 

(2) the H2 surface density, derived from CO (1-0) emission, and 
the total (H2 + HI) gas surface density are better correlated 
with the cold dust surface density than is HI alone. Thus, glob- 
ally in these galaxies, the cold dust emission is likely associated 
both with the molecular and atomic phases. 



) me r in. sunace Dngnmess increases as me mira power 
of the ratio S(60/zm)/S(100 (im). This means that for those 
objects the major part of FIR emission is thermal. The FIR 
surface brightness also shows a tight and linear correlation with 
the H2 and dust surface densities while a weaker one with the 
HI surface density. This suggests that a large part of the far- 
infrared emission of these galaxies originates in the molecular 
medium. This conclusion is stronger than what can be derived 
from the comparison of FIR luminosities and H2 masses since 
surface densities are not distance-dependent. 
(4) the mean value of the gast-to-dust ratio, (Mh 2 + Afai)/Afd 
depends on the dust model and is on average 230. Since this 
value is close to the solar neighborhood value, the estimation 
of the gas masses, (Mh 2 + Mm), is not wrong by a large fac- 
tor. It can be argued, therefore, that the computation of H2 
masses from CO (1-0) integrated intensities using the Milky 
Way empirical conversion factor is reliable. The dependence of 
gas-to-dust ratio on the FIR surface density confirms previ- 
ous suggestions of a proportionality between dust-to-gas ratio 
and heavy element abundance (Issa, McLaren and Wolfendale, 
1990; Sodroski et al., 1994). 
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Table 1 



optical and FIR properties of the sample 







IRAS position 




TV* J. 

Distance 






optical photometry 




b IK data 




name 




R.A. 
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D 


D 25 


m B 
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morp. type 
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h 


(°) 
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10 T \ 
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4.25 
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IC1553 


00 


30 10.8 


-25 


53 


02 


37.6 


70 





13.5 


2820 


SAO 


2.70 


4.81 
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11.75 


0100-22 


01 


00 22.0 
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13.6 
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4.42 


7.62 


2.88 


E478-G6 


02 


06 59.9 
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4.05 
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12300 




2.20 
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4.77 


6.17 
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11.0 
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0.37 
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Table 2 

CO and HI data 





HI data 




CO data 




name 


S 


vhi 


AV m 




vco 


AV co 




(Jy km/s) 


(km/s) 


(km/s) 


(K km/s) 


(km/s) 


(km/s) 


N142 


1.7± 0.2 


8081± 15 


250± 40 


6.2± 1.2 


8030± 40 


406± 80 


IC1553 


3.5± 0.2 


2930± 7 


151± 12 


<1 






E473-G27 


0.7± 0.2 


18338± 30 


280± 60 


3.0± 0.8 


18378± 50 


325± 80 


N232(l) 


m 2 i=16.8± 0.3 


6675± 10 


293± 11 


11.5 






0100-22 
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E475-G16 


4.6± 0.3 


7058± 12 


345± 25 


3.6± 0.9 


7018± 40 


335± 50 


N578 


m 2 i=12.83± 0.12 


1630± 4 


264± 5 


(0,0) I<0.7 














(0",45") I<1 














(0",-45") 1=1. 5± 0.5 


1665± 6 


36± 11 










(45", 0") 1=3. 3± 0.8 


1592± 9 


66± 17 


N808 


<0.5 






not observed 






E478-G6(2) 


7.5± 0.3 


5336± 15 


280± 30 


5.4± 1.8 


5352± 50 


347± 100 


N922(3) 


15.7±0.4 


3090± 3 


182± 5 


(0,0) I<1.3 






E479-G2 


2.1± 0.2 


10460± 17 


280± 40 


not observed 






N1187(4) 


m 2 i=13.03± 0.13 


1396± 5 


266± 6 


(0,0) 1=10. 9± 1.3 


1390± 5 


101± 20 










(32",-32") I=5±l 


1310± 9 


50± 20 










(64" ,-64") I<1.5 














(-32", 32") I=4.0± 0.8 


1490± 5 


50± 10 










(_64", 64") I=7±2 


1540± 16 


126±20 


N1255(5) 


20. 5± 0.3 


1685±2 


162±3 


(0,0) I=2.3±0.9 


1670±13 


70± 20 










(23", 39") 1=2. 1± 0.8 


1940± 15 


75± 20 










(-23" ,-39") I<1.2 














(-39", 23") I=4.3±1.2 


1629±15 


128±20 










(39",-23") I<1.2 






E481-G23 


7.0± 0.2 


1455± 2 


120± 4 


1.9± 0.4 


1432± 15 


93± 22 


N1385(6) 


m 2 i=13.76± 0.08 


1494±4 


181±4 


(0, 0)1=9. 0± 1.3 


1460±8 


70±15 










(43", 12") I<2.1 














(-43",-12") I=4.4± 1.7 


1498± 25 


128± 60 










(-12" ,43") I<2.5 














(12",-43") 1=3. 2± 1.1 


1561± 11 


60± 22 










(-23", 87") I<2 














(23",-87") 1=2. 9± 0.7 


1270± 6 


43± 11 


N1415 


m 2 i=15.11± 0.21 


1585± 7 


322± 5 


(0,0) 1=9. 1± 1.7 


1579± 40 


230± 60 










(32",-32") I<2 














(-32", 32") I<2 






E549-G23 


1.8± 0.2 


4226± 10 


170± 20 


2.7± 0.6 


4185± 20 


200± 40 


E483-G12 


1.5± 0.2 


4290± 30 


500± 50 


1.9± 0.7 


4247± 94 


485± 160 


N1591 


2.6± 0.3 


4166± 10 


200± 30 


3.5± 0.7 


4138± 25 


230± 30 


E484-G36 


<1 






3.8± 0.9 


5058± 17 


150± 40 


E485-G3/4 


4.0± 0.3 


4414± 7 


213± 17 


5.5± 0.9 


4564± 30 


280± 40 


IC2520 


m 2 i=15.4± 0.3 


1238± 9 


169± 5 


not observed 






MK727 


m 2 i=16.57± 0.21 


7641± 7 


177± 11 
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IC3581 


m 2 i=16.47± 0.21 


6921± 6 
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lZUlo± 4U 
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Z.Z±L U.O 
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Notes : 



- The HI data are from our observations with the Nangay radiotelescope, except when m2i is given in column 2. 
In this case the data come from the LEDA database. 

The HI masses were computed from m2i using : 
log(M ffl ) = 12.336 - 0.4 x m 21 + 2 log(D/Mpc). 

- * marks the 20% linewidths 

- ** : linewidth from a gaussian fit 

- other HI linewidths are 50% widths from LEDA. 

For our CO and HI data, the widths have been computed from the line moments. 

(1) CO data from Mirabel et al. (1990). 

(2) 2 galaxies in the beam 

(3) 5-points map at 45" sampling, I(CO) < 5.5 

(4) average value of 5 points, I(CO) = 4.1 ± 0.6. 

The mass in the Table has been extrapolated because of the large galaxy angular diameter. 

(5) average value of 5 points. I(CO) = 3.0 ± 0.5 

(6) average value of 7 points. I(CO) = 4.5 ± 0.8 

(7) data from Casoli et al. (1995, in prep.), NRAO 12m observations 

(8) average value of 3 points. I(CO) = 4.5 ± 0.8 

(9) HI and CO data from Dupraz et al. (1990) 

(10) average value of the 5 central points. I(CO) = 1.6 ± 0.4 



